Diffusion phenomena for partially dissipative hyperbolic systems 



Jens Wirth 



Abstract. 

In this note we provide some precise estimates explaining the diffusive structure of partially 
dissipative systems with time-dependent coefficients satisfying a uniform Kalman rank condition. 
Precisely, we show that under certain (natural) conditions solutions to a partially dissipative hyper- 
bolic system are asymptotically equivalent to solutions of a corresponding parabolic equation. 

The approach is based on an elliptic WKB analysis for small frequencies in combination with 
exponential stability for large frequencies due to results of Beauchard and Zuazua and arguments of 
perturbation theory. 

§1. Introduction 

The classical diffusion phenomenon observed by Hsiao and Liu [2] and later Nishihara [5] pro- 
vides an asymptotic equivalence for solutions to damped wave equations (or porous media equations) 
and corresponding solutions to the heat equation. To be precise, for any solution to the Cauchy 
problem for the damped wave equation 

(1) utt — + ut = 0, u(0, -)=ui, u t (0,-) = u 2 
on K™ we find a corresponding solution to the heat equation 

(2) w t = Aw, w(0, •) = Wq = u\ + U2, 
such that their difference satisfies 

(3) ||«(t, •)-«;(*, Oils ^^(Nlla+Klla)- 

On the other hand, both solutions do not decay in general. By this we mean, that the estimates 

(4) Kt,0l|2<C(||«i|| 3 +||« 2 || a ), IKV)ll2<C|K||2 

are both sharp. The diffusion phenomenon makes the diffusive nature of dissipation apparent. See 
also [4] and [7], [5] for a detailed study of estimates in this spirit. 

Our aim is to show that this statement has a natural counterpart in the language of differential 
hyperbolic systems with time-dependent coefficients. Our main tools will be a block-diagonalisation 
of the (full) symbol for small frequencies combined with an exponential decay rate for large fre- 
quencies deduced from methods of [5] and pQ. 
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§2. Main results 

Before stating the main results, we will introduce the classes of coefficient functions we are 
going to use later on. For a parameter r G K we denote by T{r} the class of C°°-functions on 
a : IR+ — > C satisfying the symbolic estimates 

(5) \d?a(t)\ <C k (l + t)- r - k . 

We consider a Cauchy problem of the form 

En 
k=i A k (t)B Xk U + iB(t)U, (7(0, •) = U 

for an unknown vector-valued function U : R+ x R" — > C d and with coefficient matrices 

(7) A k (t),B(t) eT{0}®C dxd 

satisfying the following assumptions: 

(Bl): the matrices A k {t) are self-adjoint; 

(B2): the matrix B{t) is non-negative and £ spec B(t) is a simple eigenvalue uniformly 
separated from the remaining part of the spectrum, 

(8) dist(0,specB(t) \{0}) > 7 > 0; 
and 

(B3): the matrices B{t) and A(t,^) = Yjk=i A k (t)£ k satisfy 

(9) -M^^rl^llBm^t^M^cin 2 , t>t Q 

for any choice of numbers ao, . . . , ad—i > and suitable constants c and to depending on 
them. 

Assumption (Bl) guarantees that the system ([5]) is hyperbolic. The non-negativity of B{t) implies 
dissipativity. The non-triviality of the null space of B(t) in (B2) is not the crucial assumption as this 
can always be achieved by multiplying U (t, x) with a scalar function and the main result remains 
true modulo the obvious modifications. Crucial point for our result is that the lowest eigenvalue is 
uniformly separated from the others. 

The last assumption is used to treat non-zero frequencies. We refer to (B3) as uniform Kalman 
rank condition. It allows to construct a Lyapunov functional in order to prove that the Fourier 
transform of solutions to (|6]) with respect to spacial variables satisfies 

(10) ||t/(t,-)l| 2 <Ce- ct K] 2 ||C/o|| 2 , [£] = |£|/<£)^min{|ai}, 

with suitable non- negative constants C, c > 0. The proof for this fact is essentially the same one 
given in [1] with the obvious modification that the Lyapunov functional has coefficients in 7~{0} 
such that heir argument works only if t is large. 

If assumptions (B1)-(B3) arc satisfied then we can associate to ^ a parabolic equation 

(11) d t w = V • a(t)Vw + pit) ■ Vw + j(t)w, w(t ,-) = wo 

with coefficients a(t) G T{0} <g> C dxd , /3(t) e T{1} ® C d and j(t) e T{1}, a smoothing linear 
operator mapping W : Uq <— > wq and a vector K(t, D) of second order differential operators with 
coefficients in T{0} such that the following theorem holds true. 
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Theorem 1. For any solution U(t,x) to © to initial data Uq G L 2 (M. n ;C d ) the corresponding 
solution w(t, x) to (jlip with wq = WUq satisfies 

(12) ||l7(t,-)-«'(t,D)«;(t,.)||2<Gt- 1/a log(t)||l7o||2 ) t>t Q . 

The coefficients of the parabolic problem (fTTj) as well as of the operators W and K(t, D) will 
be made precise in the following section. They arise naturally within a decoupling procedure of the 
system. 

§3. Proof 

The proof is divided into several parts. First we will asymptotically block-diagonalise the full 
symbol of the problem (J6]) for small frequencies within the so-called elliptic zone Z c \\(e) = {(£,£) : 
t > e _1 , |£| < e} for sufficiently small e. Similar to [ID] , this will allow us to obtain asymptotic 
formulas for the fundamental solution for bounded t\£\ < 5 within this zone. In a second step we 
will use main terms of this asymptotic representation to deduce the coefficients of (|11|) . In a final 
step we will prove the desired estimate of Theorem [1] chosing k and S sufficiently large. For full 
details and further information we refer to [9[ Chapter 5] . 

3.1. Block-diagonalising A(t, £) + B(t) 

The consideration follows partly [3] and |TTj . In the sequel we will use the notation 

(13) r{m} = L( t ,^ = y2 Pa (t)C-P a (t)€T{m-\a\}\ 

* — '\a\<m 

for polynomials of degree m with coefficients in the 7~-classes. We will use the same notation for 
the matrix- valued case. Clearly V{m + 1} C V{m} such that these classes constitute a hierarchy. 
Condition (JT]) guarantees that A(t, £) G ^{l}, while B(t) G V{Q}. Furthermore, due to © we also 
know that there exists an invertible matrix M(t) G V{0} such that 

(14) M-\t)B(t)M(t) = b-diag^^i) (0, B(t)) = V{t) 

is block-diagonal. Denoting V(°'(t,£) = M _1 (t)[/(t, £), U being the partial Fourier transform of U 
with respect to spatial variables, we obtain for the ordinary differential equation 

(15) D t V<°> = (V(t) + i?!(t,0)^ (0) 

with i?i(i,0 = AI- 1 (t)A(t,£)M(t) + (ptM^^M^) G V{1} of lower order in the ^-hierarchy. 
By induction we show the following lemma. 

Lemma 2. Let k G N, k > 1. Then there exists a constant Ck and matrices Nk(t,£) G V{0}, 
block- diagonal F k (t,^) G and Rk+i(t,^) G V{k + 1} such that 

(16) (D t -V(t)- R 1 (t, 0) N k (t, = N k (t, £) (D t - V(t) - F k (t, - R k+1 (t, £)) 



holds true within Z e \\(ck)- The matrix Nk(t,£,) is uniformly invertible within this zone. 
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Proof. We explain the k = 1 step in detail. We construct a matrix N^(t,^) 6 V{1} such 

that 

(17) (D t -V(t)-R 1 (t,O)0. + N^(t,O) 

- (I + NV>(t,£))(D t ~ V{t,0 - F 1 (t,0) 6 ^{2} 

holds true for some diagonal matrix Fi(t, £) G ^{l}. Collecting all terms not belonging to the right 
class yields again conditions for the matrices N^\t,^) and F\(t, £). Indeed, 

(18) [D(t),NW(t,0} = -R 1 (t,0 + F 1 (t,0 
must be satisfied and therefore it is reasonable to set 

(19) F 1 (t,0=b-diag {lid _ 1) R 1 (t,0 

and to denote by NW(t,£) = ( ° (< ?) rll . 1 (*^) T J the matrix with entries 

(20) n 1 , j (t,0=J o °°e- s ^R 1 , j (t,Ods, R.-F^^^ R ^ , 

solving Sylvester's equation (TPS]) . This implies that iV^ 1 )^, £) € Recursively, we will 

then construct matrices N( k '(t,£) £ V{k} and F"(t,() S ^{^l block-diagonal, such that for 

(21) N K (t, o = i+ J2 K k=1 Nik) (*. o = EL ^ (fc) (*> o. 

the estimate 

(22) B K (t,0 = {D t -V(t) - J2 1 (* > e))JVjc(*,0 

- Nx(t, 0(D t - 2>(t) - 0) G ^{/v + 1} 

is valid. As we just did this for = 1, it remains to do the recursion k H >• fc + 1. Assume 
The requirement to be met is that 

(23) B k+1 (t, - B*(t, = - [X>(t), A( fc+1 ) (t, f)] + ^ fc+1 > (t, mod + 2}, 

which yields again a Sylvester equation and its solution is given in analogy to the above case. It is 
evident that the construction implies (t, £), jV( fc+1 ) (t, £) £ V{k+1} together with B k+1 (t, f) S 

7>{fc + 2}. The matrices jVk(i, £) £ P{0} are invertible with inverse N^ty, f) £ V{0} if we restrict 
our consideration to a sufficiently small elliptic zone Z e \\(c k ). Q.E.D. 

3.2. Asymptotic integration 

We consider the transformed problem in I^ (fc )(i,£) = N k (t, (t, £), 

(24) D t yW(t,o = (v(t)+F k (t,t) + R k+1 (t,o)vW(t,o, 

and reformulate this as integral equation for its fundamental solution. Denoting it by £ k (t, s, £), we 
know that it solves the above equation to initial data £ k (s, s, £) = I £ C dxd . 
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Let &k(t, s, be the fundamental solution to the block-diagonalised system D t ~V{t)—Fk{t 1 £). 
Then 

(25) e k (t,s,o = (^ k{t s '° e fc (M,o)' H^(^ s ^)ll<c-^ 

holds true for t > s uniformly within Z e \\(e) for e < c k sufficiently small. Furthermore, £k(t,s,£) 
satisfies 

(26) £ k (t,s,0 = Q k (t,s,Z) + J O k (t,9,Z)R k+1 (0,O£k(O, S ,OdO. 
We solve this Volterra integral equation using the Neumann series 

(27) f fc (t,a,0 = Qfc(* ) *,0 + X)Zi/ e*(*.*i.Ofl*+i(*i»0 ••• 

■ ■ • y fc (t/_i, t/, £)Rk+i(ti, Odtt ■■■dt 1 . 

This series converges based on the estimates for the remainder term iJfc + i(t, £) 6 T'{k+ 1}. It can 
be estimated by 

(28) \Mt,s,0\\ < exp (J \\R k+l (e,0\\d9^ . 

The series converges uniformly within Z e \\(c k )C\{t\£,\( k+1 ^ 2 < 8} for any constant S and ||£fc(f , s, £) — 
Q k (t, s, £)|| — > as c k — > for fixed (5 > as soon as we choose k > 1. 

We can do slightly better. Let VF(s) = lim^oo ej£ k (t, s, 0). Then the following lemma holds 
true. The full proof can be found in [5] . 

Lemma 3. The fundamental solution £ k (t,s,£), k sufficiently large, satisfies the estimate 

(29) \\£ k (t,s,0-Z k {t,s,£)eiW 2 (s)\\ <C k (l+t)~i 
uniformly in |£| < e k . 

3.3. The parabolic reference problem 

The parabolic reference problem is defined in terms of the upper left corner entry of the matrix 
F k (t,£) modulo terms from P{3}, 

(30) f[ k \t,0 = i]T rf a id (tMj + + mod ^{3} 

with Q!ij(t) G T{0}, € T{0} and 7(i) G T{1}. Now the uniform Kalman rank condition (B3) 
implies a structural property. It follows that (3i(t) has to decay, 0i(t) G 7~{1}, and also that the 
quadratic matrix {ctij(t))ij is positive definite modulo 7~{1}. The latter is a direct consequence of 
the estimate (flT)|) . 

This defines the coefficients a(t), f3(t) and ^(t) in (jTTJ) . The corner entry E2(t, s, £) of 62(^7 s, £) 
is the fundamental solution to this Cauchy problem. Furthermore, we define K(t, £) = M(t)N 2 (t, £)ei, 
where d = (1, 0, . . . , 0) T G C d . By definition we have K(t, £) G V{2}. 
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3.4. The final estimate 

We define wo = WUq in such a way that we cancel the main term of the solution within 
-Zeii(cfc) n {t\£\ < 5}, i.e., we define 

(31) = W 2 (t 0l ON2 1 (toA)M- 1 (t )£(t 0l O,Ox(OUo 

with x(0 £ Co°(M"), x(0 = 1 ncar £ = and suppx C B C2 (0). Then the estimates of Lemma[3]im- 
ply the following statement. The logarithmic term is caused by comparing B,k(t, s, £) with Ea(t, s, £). 

Lemma 4. TTie solution w(t, x) to (| 1 1|) satisfies 

(32) ||% •) - A'(i, •)©(*, -)||a < C"(l + t)- 1 / 2 log(c + 1). 

§4. Concluding remarks 

We used that ker B(t) is one-dimensional. The result can be generalised to situations where a 
(stable) higher-dimensional null-space appears and further diagonalisability conditions on appearing 
lower order terms in the diagonalisation hierarchy are satisfied. This corresponds to the block- 
diagonalisation in [3] and again yields a parabolic reference system for these diffusive modes. 

If assumption (B3) is violated, a variety of other asymptotic scenarios may appear. If n = 1, 
the results of [I] hint to a decomposition of solutions into traveling waves, parabolic type modes 
and exponentially decaying modes. 
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